IPPP-01-25 
DCTP-01-50 
bep-ph/0205l20 
May 2002 



SEMI NUMERICAL RESUMMATION OF EVENT SHAPES AND JET RATES 

GIULIA ZANDERIGHI 

Institute of Particle Physics Phenomenology, University of Durham, Durham DH1 3LE, England 



For many observables, the most difficult part of a single logarithmic resummation is the ana- 
lytical treatment of the observable's dependence on multiple emissions. We present a general 
numerical method, which allows the resummation specifically of these single logarithms. Some 
first applications and new results for the thrust major, the oblateness and the two-jet rate in 
the Durham algorithm are also presented. 

1 Introduction 

The formation of jets in hadronic collisions is surely one of the most striking phenomena in 
high energy collisions. Their study allows a precise measurement of a s , is important for any 
search for new physics and is useful to test our understanding of strong interaction dynamics. 
Generally jet cross sections are quite large but it might be difficult to get an insight into details 
of the interaction mechanics. One way is provided by the study of event shapes and jet rates. 

2 Resummation of event shape distributions 

Event shape variables can be studied at variety of levels, either through relatively inclusive 
properties, such as their mean value, or more exclusively by examining integrated (or differential) 
distributions. The integrated distribution accumulates contributions both from real and 

virtual emissions. In spite of the fact that both contributions are separately divergent, infrared 
safety of the observable ensures that their sum is finite at any order in perturbation theory. In the 
(more inclusive) v ~ 1 region the observable is sensitive only to the total energy and momentum 
flow, while in the (less inclusive) region v — * the mismatch between real and virtual emissions 
gets emphasized, so that divergences do still cancel but large logarithmic enhanced contributions 
are produced, which need to be resummed at any order in perturbation theory. 



"From here onwards, capital letters denote the observable, while lower case the numerical value. 



2. 1 Classification of single logarithmic terms 



In spite of the fact that resummation programmes are quite cumbersome it is possible to classify 
systematically in a simple way all leading (DL) and subleading (SL) contributions. DL contri- 
butions are due to emissions which are both soft and collinear, while SL corrections are found 
to be due to 

• soft, large angle emissions; 

• hard, collinear radiation; 

• effects due to the running of a s ; 

• multiple emission effects, i.e. taking into account properly how all emissions add up 
together and interfere to contribute to the value of the observable. 

As far as the first three items are concerned, to obtain a resummed prediction at SL level it 
is enough to understand how a single emission affects the observable and to exponentiate this 
result in a naive way. Therefore the last point turns out to be the only non trivial task in any 
resummation programme and it is exactly this contribution which we address here. 

2. 2 Standard treatment of multiple emission effects 

To understand how multiple emission effects are treated usually let us consider for instance the 
single jet broadening B (in the following ku denote the secondary transverse momenta) 

2BQ = ^2k t i, with k ti = \k ti \ . (1) 

i 

Resummation is achieved factorizing matrix element and the phase space (using an independent 
emission picture □'□) |] and the definition of the observable via a Mellin transform 

Q{2BQ - ku) = jf -^~/ VBQ z- vku • (2) 

But, actually, at SL level one needs to take into account the fact that due to all secondary 
emissions the hard parton will take a recoil pt- Therefore the 0-function needs to be replaced 
by 

®{2BQ - hi) => J d 2 Pt S 2 (p t + £ k ti ) Q(2BQ - £ Aft - Pt) ■ (3) 

i i i 

To factorize the 5-function one needs now an additional two-dimensional Fourier transform. In 
the case of the thrust minora, which is defined as the radiation out of a plane, to fix kinematically 
the plane one needs an additional five-dimensional integral, so that generally an analytical 
treatment can become very involved. Furthermore, in certain cases, as for the thrust major, the 
oblateness and the 3-jet resolution, an analytical treatment seems to be unfeasible. 

3 Numerical treatment of multiple emission effects 

Since multiple emission effects have the same origin it is natural to look for a general approach 
to the problem. The idea is to relate multiple emission effects of the observable V under study 
to a 'simple' reference variable V s which has the same double logarithmic structure, i.e. the 
same dependence on one soft-collinear emission, but which factorizes in a trivial way. 

'Actually, the single jet broadening, taken here for illustration, is a non-global observable, so that an indepen- 
dent emission picture mistreats some of the single logstj. 



For an observable given by the sum of n contributions, i.e. sensitive to all emissions, the 
natural choice of the simple variable becomes 

V = Y,v(ki) V s = max[v(ki),...,v(k n )} . (4) 

i 

By construction V s is sensitive only to the 'largest' emission and this definition also ensures that 
V s factorizes in a simple way (i.e. needs no Mellin transforms) 

e(V s -m a x[v(k 1 ),...,v(k n )])=l[e(V s -v(k i )). (5) 

i 

As a consequence the resummation of V s becomes straightforward. For instance for the case of 
e + e~ — > 2jets the resummed distribution becomes (the case e + e~ — > 3jets is also trivial) 



X s (v s )=e- H ^°> K s (v s ) = C F J / Qe : 3 /4 dy @(v(k) - v.) ■ (6) 

* fc t 

The aim is now to exploit the 'simple', known distribution S s to compute the real observable E. 
The two differential distributions D = dH/dhxv are related by a simple convolution 

D(v) f dv. 



f ^l Ds (v s )P(v\v s ) , (7) 
J v s 



where P(y\v s ) denotes the conditional probability of having a value of v for V given a value of 
v s for V s . The two distributions have the same double logarithmic structure so that v ~ v s and 
one can expand D s (v s ) ~ e - ^ 8 ^ 8 -' around the value of v s 

D s (v s )= SL D s (v)e- R ' ln ^ , (8) 

where R' = —dR/dlnv is a single logarithmic function. Combing eq. @ and @ one obtains 

D{v) = SL D s {v)HR') , (9) 

with 

F{R') = { — e - R ' lnx p(x,R') , x = v/v s , p[x,R')=vP(v/v s ). (10) 
j x 

At SL level the same relation holds then for the integrated distribution 

E(t;) = ^ s (v)T(R') . (11) 

This expression shows explicitly that all non trivial multiple effects can be factored out and 
embodied in T. The aim is then to compute T i.e. p(x, R') in a general way. 

3. 1 Numerical implementation: the procedure to get J- 

To compute the function T numerically one writes a Monte Carlo! which starts from an arbitrary 
Born configuration with an arbitrary number of hard emitting legs and 

generates the first 'largest' emission with v{k\) = v s (i.e. one starts fixing V s = v s ); 

1 generates an additional soft-collinear (SC) emission (according to the phase space R') 
which satisfies v(k{) < t>(/cj_i); 

iT if v(ki) < e <C 1, where e defines the precision of the procedure, the emission (and all 
subsequent ones) is too small to affect the value of the observable, therefore the procedure 
is stopped and the momenta generated are passed to a standard routine which computes 
the value of V; otherwise one sets i — > i + 1 and goes back to step /. 



4 Recollection of known results and new results 



As a first check of the procedure we compared some Monte Carlo results with results from the 
literature for the function T for the thrust u and the total and wide broadening □ and for the 
thrust minor §. Comparisons show that our numerical results turn out to be interchangeable 
with analytical ones til. 

We then applied the procedure to obtain some new predictions for thrust major and the 
oblateness, which have so far defied analytical resummation and to the two-jet rate in the 
Durham algorithm, for which only a subset of the single logs had up to now been calculated B0. 

We first used the Monte Carlo procedure to compute the function T for these variables, 
combining this with £ s we get the first full resummed prediction at SL level which we then 
match with second order exact results. 




5 Conclusions 

The study of event shapes and jet rates has proved to be a powerful tool to test QCD and to 
measure a s . Many variables have been studied in the last years, but never in a general way. 
This work is a first step towards a complete numerical resummation of any n-jet observable in 
any hard process. What has been done so far is to implement multiple emission effects generally, 
to test the procedure and to obtain some new predictions. What remains to be done is to make 
the procedure completely automatic, i.e. userfriendly and to exploit the method to obtain a 
bunch of new predictions especially in hadron-hadron collisions, which will be the main object 
of most experimental analyses in the near future. 
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